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We study the effect of dense quarks in a SU (N) matrix model of deconfinement. For three or more 
colors, the quark contribution to the loop potential is complex. After adding the charge conjugate 
loop, the measure of the matrix integral is real, but not positive definite. In a matrix model, quarks 
act like a background Z(N) field; at nonzero density, the background field also has an imaginary 
part, proportional to the imaginary part of the loop. Consequently, while the expectation values 
of the loop and its complex conjugate are both real, they are not equal. These results suggest a 
possible approach to the fermion sign problem in lattice QCD. 

PACS numbers: 



At nonzero temperature, numerical simulations in lat- 
tice QCD have provided fundamental insight into the 
transition from a hadronic, to a deconfined, chirally sym- 
metric plasma At nonzero quark density, however, 
at present simulations are stymied by the "fermion sign 

problem" H El a i i S i El M El El El 13 EE 

Il7j . Even in the limit of high temperature, and small 
chemical potential, only approximate methods can be 

used mm Ei mm. 

In this paper we consider deconfinement in a mean field 
approximation to a model of thermal Wilson lines I2J, 
p| .which is a matrix model [II El El El El El III 
E1E1E1E1E3- In Sec. we discuss general features of 
SU (N) matrix models at nonzero quark density j2^] • In 
scc.|n| this is briefly contrasted with the (trivial) case of 
a U(l) model Numerical results for three colors are 
presented in Sec. IIIII In Sec. IIVI we conclude with some 
remarks about some methods which might be of use for 
dense quarks in lattice QCD. 



I. SU(N) MATRIX MODEL 

In a gauge theory at nonzero temperature, a 
basic quantity is the thermal Wilson line, L = 
Vexp(ig J A dr), where g is the gauge coupling, A 
is the timelike component of the vector potential, and 
the integral over the imaginary time, r, runs from 
to l/T, where T is the temperature [24[. An effec- 
tive theory of thermal Wilson lines, interacting with 
static magnetic fields, can be constructed, and is valid 
in describing correlations ove r sp atial distances ^> l/T 

IH El El IS El El El El EI E3 

Over large distances, we use a mean field approxima- 
tion to this effective theory. This gives an integral over 
a single Wilson line, L, with the partition function that 



of a matrix model: 

Z = J dL exp (-(A 2 - 1) (V fl i(L) + V, fc (L))) ; (1) 

L is an SU(N) matrix, satisfying L^L = 1 and det L = 1. 
Under gauge transformations Q, it transforms as L — » 
Q^L£7, so that gauge invariant quantities are formed by 
taking traces of L. These are Polyakov loops. In the 
matrix model, the effects of gluons and quarks are rep- 
resented by potentials, V g i(L) and Vqfc(L), which are 
(gauge invariant) functions of the Wilson line. The ef- 
fects of fluctuations, which are not included in the ma- 
trix model, can also be included in a systematic fashion 

[HEi 

The pure glue theory is invariant under a global sym- 
metry of Z(N) , and so this must be a symmetry of the 
gluon loop potential, V g i(L). The simplest form for the 
gluon loop potential is a type of mass term, 

V a i(L) = -m 2 £ w £ N , ^v = ltrL, (2) 

where £jv, and ijf = £* N , are the Polyakov loops in the 
fundamental, and anti-fundamental, representations. Up 
to a constant, this gluon potential is proportional to the 
Polyakov loop in the adjoint representation. 

In general, the gluon potential is a sum over all loops 
in Z(N) neutral representations |3^|; this can be written 
as a power series in terms like (|^tv| 2 ) 2 , etc. These terms 
are invariant under a larger global symmetry of U(l). 
The first term which is invariant under Z(N) , but not 
U(l), is 

(i N ) N + (e w f. (3) 

Another such term is 

H(eN) N -(e w ) N ) , (4) 

where the factor of i is added to ensure that in all, the 
term is real. 
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While (JSJ certainly appears in the gluon loop poten- 
tial, terms such as Q should not arise in effective the- 
ories of relevance to QCD. Gluons are invariant under 
the discrete symmetry of charge conjugation, C, under 
which — ► — A* (takin g A u , = A^t a , and Hermitean 
generators t a for SU(N)) g^. Under C, the Wilson line 
transforms into its complex conjugate, L — > L*, so that 
© is even under C, and Q}, odd. 

Quarks in the fundamental representation of SU(N) 
are not invariant under the global Z{N) symmetry. Thus 
quarks tend to induce a background Z(N) magnetic field, 
which we characterize by a parameter h. The simplest 
contribution to the quark loop potential is then |2y, 1271 

V,*(L) = ~ (e" ^ + c-"%) . (5) 

At finite N, h affects the deconfining transition 
in the standard manner of a background magnetic field 
[2(jl l30l |. If the deconfining transition is of first order 
in the absence of quarks, then their presence tends to 
weaken the transition. Eventually, it disappears at a crit- 
ical end-point, for some value of h; above this value, there 
is no phase transition, just a smooth crossover. If the de- 
confining transition is of second order in the absence of 
quarks, then any background field, h ^ 0, washes out the 
transition. 

At infinite N, if one is away from the Gross-Witten 
point, then the behavior is like that at finite N. Precisely 
at the Gross-Witten point (3^, , correlation lengths 
diverge at the transition. Then like a second order tran- 
sition, any background field changes the order: from first 
order at the Gross-Witten point, into one of third order 
when h^O [HUH. 

We have also added a parameter, fi, to represent the 
quark chemical potential; fi should be understood as the 
true quark chemical potential, divided by temperature. 
The quark chemical potential is associated with a con- 
served charge for the global [7(1) symmetry of baryon 
number. This dictates that the chemical potential enters 
in the above form, like the imaginary component of a 
[7(1) gauge field 0. 

Under charge conjugation, the Wilson line transforms 
into its complex conjugate, and the chemical potential 
changes sign: 

C : L -> L* , n -y -n . (6) 

The term in JSJ is invariant under C, as should be all 
terms in the quark loop potential. 

Implicitly, we have integrated out the quarks to ob- 
tain the loop potential in (jSJ). For example, if one com- 
putes the quark determinant in a background gauge field, 
~ tr log(^>+ m q k), one will obtain a term such as JSJ: 
see, for example, the calculations of Langfeld and Shin 
[l2| and of Schnitzer [3(j . Other discussions of l oop po- 
tentials with quarks include those of [2^, |2^, |3(], |3l| ; at 
nonzero quark density, see [27l l3l| . These calculations 



show that at a temperature T, the background field of 
massive quarks behaves as h ~ exp(— m q k/T), reaching 
some finite value as the quark mass vanishes. As with the 
gluon loop potential in J2J, there are many other terms 
besides that of in V q k- These involve all possible 
traces of 

c' 1 L and e^ L* . (7) 

in such combinations which are invariant under charge 
conjugation, ©. These two matrices represent, respec- 
tively, the propagation of a particle forward in imaginary 
time, and an anti-particle backward in time. Of course, 
charge conjugation symmetry is violated by a given value 
of /i =/= 0: C just implies that a Fermi sea of quarks be- 
haves similarly to one of anti-quarks (neglecting electro- 
weak interactions). 

The quark contribution to the loop potential equals 

V q k(L) = — h (cosh(/x) Re £n + i sinh(/u) Im £n) , (8) 

where Re and Im denote the real and imaginary parts, re- 
spectively. At zero chemical potential, quarks generate a 
real background Z(N) field for the real component of the 
loop, ~ Re^Ar. When the chemical potential is nonzero, 
however, the background Z(N) field not only contains 
a piece proportional to the imaginary part of the loop, 
~ Im £jsf , but with a coefficient which is itself imaginary. 

The case of two colors is special. For two colors, loops 
in any representation are real, and for any /i, the back- 
ground field generated by quarks is always real. For three 
or more colors, however, loops have imaginary parts, and 
the potential generated by quarks is manifestly complex, 
(jHJ). This is how the fermion sign problem appears in a 
matrix model. 

In this case, though, it is easy to reduce the sign prob- 
lem, which appears to be one of complex phases, to one 
in which the phases are always real. If a given matrix, 
L, contributes to the partition function, then so does its 
charge conjugate, L*. Adding the contributions of L and 
L* together, we can rewrite the partition function in a 
form which is manifestly real, 

Z = j dLe- v(L) cos(^Im£Ar) , (9) 

where 

V(L) = {N 2 - 1) (V g t(L) - h cosh( A1 ) Re l N ) , (10) 

h = (N 2 -l)h sinh(^) . (11) 

The potential V(L) is even under charge conjugation of 

the gluons, while h Im£jv is odd. We can use this to 
write the expectation value of the fundamental loop as 

(£ N ) = A JdL e"^ L ) x 
(^cos(hlm£ N ) Rc£ N — sin(hlm£ N ) Im^^ ,(12) 
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while that of the charge conjugate loop is 



1 



dL e"^ 



Nl z 

(cos(h lm£pf) Re^Ar + sin(/i Imfjv) Im^jvJ .(13) 

Because dense quarks induce an imaginary background 
field for the imaginary part of £n, the expectation values 
of £jv and lj? are not equal to one another, although they 
are both real. 

Physically, this is natural. A loop is proportional to 
the (trace of the) wave function of a quark; the complex 
conjugate loop, to that of an anti-quark. A Fermi sea 
represents a net excess of quarks over anti-quarks, so at 
/i ^ 0, quarks and anti-quarks propagate differently. In a 
matrix model, this manifests itself as unequal expectation 
values for £n and £- N -. 

Karsch and Wyld performed numerical simulations for 
a model of SU(3) matrices, living on sites of a three- 
dimensional lattice, at nonzero density Q. Our matrix 
model represents a mean field approximation to their the- 
ory. They were the first to observe that the expectation 
values of £n and £jj differ at nonzero density. This also 
happens for a Potts model at nonzero density 0, ^| . 

This contrasts with what would happen if the back- 
ground field which coupled to the imaginary part of L 
was real; i.e., for \x = ifi. This corresponds to a U(l) 
rotation of L, so that both expectation values are com- 
plex, and satisfy (In) = ((%-})*■ In a U(N) theory, this 
just rotates the vacuum by an angle = pt; for SU(N), 
because of the Z(N) symmetry, the vacuum structure is 
more involved. 

In Sec. IIIII we present numerical calculations of 
the expectation values of the fundamental and anti- 
fundamental loops for N = 3. Even without explicit 
calculation, however, we can understand the qualitative 
nature of the solutions. 

Consider first the limit about zero chemical potential. 
Taking the derivatives of the expectation values in 112|) 
and ifT^jl with respect to fi, 



d(£ N ) 



dfi 



d(£- 



Nl 



Of I 



= - h{N 2 -l)((lm £) 2 ) 



fj,=a 



as y, increases. 



so does 



}, 



(14) 
while 



About /i = 0, then 
(£n) decreases. 

It is also easy to understand the behavior of the ex- 
pectation values in the limit of large /i. This corresponds 
to a very strong background field, proportional to ~ £n- 
Taking the Wilson line L = exp(iw), the real part of £n 
is ~ tr lu 2 , while the imaginary part is ~ tr to 3 . For large 
background field, then, the potential is dominated by the 
real part, ~ ftexp(^i)tr w 2 ; fluctuations in the imaginary 
part are suppressed, by ~ exp(— fi/2) relative to the real 
part. Thus as /i — > oo, the expectation values of £jq and 
£jj both approach unity, (£ N ) w (£jj) — > 1. 

(The parameter \i is the quark chemical potential di- 
vided by temperature, so naively, T — > corresponds to 



/i — > oo. Remember, though, that our effective theory is 
only valid for distances 3> 1/T. We believe that the large 
/i behavior is an artifact of the model, and is not indica- 
tive what happens in the full theory at low temperature. 
See, also, Sec. IV.) 

We can thus anticipate the behavior of the expectation 
values of the loops as a function of [i. Due to the back- 
ground Z(N) field of the quarks, both loops are equal 
at fi = 0. As /i increases, at first the two expectation 
values split: one increases, while the other decreases. As 
fj, — ► oo, they come together and approach unity. For 
N = 3, this is illustrated in Sec.lrTTlbv Fig.|U 

It is customary to interpret the expectation value of 
the Polyakov loop as the "free energy" of a test quark 
|25j . At nonzero density, this implies that the expectation 
value of the fundamental loop is the "free energy" of a 
test quark, and that of the conjugate loop is the "free 
energy" of an anti-quark 

(£ N ) = exp(-F,/T) , (£ w ) = exp(-F ? /T) . (15) 

Any free energy, however, should decrease monotonically 
with /j,; because (£n) decreases about /x = 0, though, 
the "free energy" for a test quark increases with /j,. This 
quandry is resolved by recognizing that the expectation 
values of the loops are not free energies, but just the 
traces of test propagators [3 |33 . As such, they need 
not behave monotonically with [i. 



II. [/(I) MATRIX MODEL 

Before going into numerical results for SU(3), we 
briefly discuss what happens in a U(X) model, as first 
proposed by Gibbs [5|. 

For U(l), the loop is just £ = exp(i6), where 9 runs 
from — 7r to 7r. At a nonzero density /i, we take the par- 
tition function as 



Z 



d6 exp ( — (e 



*£*) 



(16) 



Like SU (N) , the fermion contribution to the loop poten- 
tial is complex at nonzero chemical potential. Summing 
over a given 9, plus its charge conjugate, which is just 
—9, the partition function becomes 



Z = 



d9 ' 



sh Wcos9 cos(/isinhOu)sin0) , (17) 



which is real. 

As for SU(N), the expectation value of a loop, and its 
charge conjugate, are unequal when /j ^ 0. However, in 
the original integral, (|16f) . we can shift the integration by 



i/j, 



(18) 



Doing so, we find that the partition function is com- 
pletely independent of fi. In terms of expectation values, 
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this implies that 



(19) 



This is an immediate consequence of the change in vari- 
ables possible in a U(l) model, (|18J) . 

For SU (N) loops, we found that both loops approach 
unity at large \i. This is not true for U(l) loops, (jT§|l : 
as /j, — > oo, ( i ) is very small, while ( £* ) is large. The 
difference arises because for J7(l), the real part of the 
loop is cos#, while the imaginary part is sin (9. At large 

the real part is ~ 1 — ff 2 /2 + . . ., while the imaginary 
part is ~ 9. At large fi, then, for U(l) the imaginary 
part of the loop dominates, instead of the real part, as 
for SU{N). 

There is a simple physical reason why, for t/(l), the 
partition function is independent of the fermion chemical 
potential With a U(l) gauge field, there is no way of 
forming baryons: the only states which are neutral under 
U(l) are trivial, having an equal number of fermions, 
exp(i9), and anti- fermions, exp(— iff). 

There is a less trivial consequence of this observa- 
tion. Consider a SU(N) gauge theory. Up to correc- 
tions ~ 1/iV, at large N, there is no difference between 
the measure for SU(N) and that for U(N). For a U(N) 
gauge theory, however, we can rotate the quark chemi- 
cal potential away. To the extent that SU(N) is close 
to U(N), then, at large N the effects of quark chemical 
potential appear only in terms which are subleading. 

Put more directly, assume that deconfinement, and chi- 
ral symmetry restoration, occurs at some temperature 
Td ~ T x when fj, — 0. Then the natural scale at which 
the quark chemical potential matters is larger than Td by 
some (fractional ) po wer of N, which can be computed in 
a matrix model 1381. 



III. N = 3 MATRIX MODEL 

In this section we present numerical results for three 
colors, where the SU(3) matrix model is just a two di- 
mensional integral. When the chemical potential \i is 
real, and the background field h is large, the integrands 
of ©, (|12p. and (|13[1 oscillate strongly. Nevertheless, we 
show that the value of these integrals are not sensitive 
to large cancellations of positive and negative contribu- 
tions, and can be computed numerically without great 
difficulty. 

For three colors, the loop potential is a function of the 
triplet and anti-triplet loops, 



1 



tr L 



4 



1 



tr V 



(20) 



We straightforwardly extend the analysis of j3j|, going 
into some detail in order to avoid confusion. Previously, 
we assumed that the expectation values of the triplet and 
anti-triplet loops are equal; now we must allow that they 
can differ. In the partition function of we introduce 



two fields, A and A*, which are the values of these loops 
for a given matrix, L: 

Z = J dhj d\* J d\ S(K - %) 5{X - 4) 

e X p(-8(V gi (A*A)+V gfe (A*,A))) . (21) 

We then exponentiate the constraints by introducing 
fields ZU* and ZU, 

Z = y e?A* J doJ* J d\ J duJ J dL exp (— 8V CO ns.) , 



Vcons. = V g i + V qk + iZJ* (A„ - %) + itj(X - 4) . (22) 

At all stationary points iZJ and ioj* are real, so we define 
lo = iZJ and oj* = iSJ*. 

Next, we define the matrix integral 



■-gw\ 



(w* ,w)= / dL exp ( 8 (w* % + w l 3 )) . (23) 



For given values of and u>, the expectation values of 
the fields are 

^!(u>„w) = i J dL % exp (8 (u £ 3 + %)) , 



dL 4 exp (8 (to 4 + %)) . (24) 



We introduce the Vandermonde potential, as a function 
of two fields, I and £*, through Legendre transformation, 



Zgw{u*,lu) = J d£* J di 

x exp(8(w»4 + lu£ - Vvdm (4, £))). (25) 
The stationary point of this integral is for 



to*(£°,e°) = 



SVvdm(t.,l) 



(£,=i° ,1=1°) 



di 



(£,=£° ,i=e°) 
This satisfies the consistency condition 
dw.{l°,iP) du(C£°) 



(26) 



(27) 



For given values of and to, we numerically computed 
the integrals in l)24|) , to obtain £° and £° . We then invert 
them, to obtain uj* and w, as a function of £° and tr. The 
Vandermonde potential then follows: 



Vvdm(£*,£) = J dt° lu(0,£°) + J dtl u.(g,i) . (28) 
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We have chosen a definite path to go from (0, 0) to (£*,£), 
but because of (|27|) . the integral is independent of the 
path chosen. 

The complete effective potential is the sum of the 
gluon, quark, and Vandermonde potentials: 

Ve// = Vgl + Vqk + Vvdm • (29) 

The Vandermonde potential, Vvdm (£* > £) , represents the 
effects of the SU(3) measure, and so is invariant un- 
der Z(3) transformations, I — » exp(27d/3)£ and — > 
exp(— 27ri/3)€». In contrast, the quark loop potential is 
not Z{3) invariant. 



To illustrate the physics, in Fig. ^ we show three ex- 
amples, with ft = 0.0 or 0.1, and \£\ = 0.2. When there 
is no background field, ft = 0, there are three degen- 
erate minima at 9 = and ±27r/3. When ft ^ and 
Jl = — 27r/3, the background field "tilts" the potential so 
that the expectation value is along the opposite direction, 
for 9 = 2tt/3. Lastly, when ft, ^ 0, and for the special 
choice of Jl = 7r/3, the background field points exactly 
in the direction opposite to the minimum at 27r/3; then 
there are two degenerate minima, for 9 — and —2tt/3. 
The potential for other values of 9 and I follows similarly. 
Also, it is clear that, as a function of Jl, the expectation 
value of 9 is discontinuous at Jl = 7r/3, jumping from to 
— 27r/3. Analytic continuation of Z to real fi is therefore 
possible only for < n/3. 
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FIG. 1: The effective potential V e ff for the SU(3) matrix 
model at imaginary /x; there is no gluon loop potential. In 
all curves, \£\ — 0.2, with 6 the phase of the loop: the solid 
curve is h = 0, the dotted curve ft = 0.1 and /i = — 2ivi/3, the 
dashed curve h = 0.1 and fi = i7r/3. 

As a check on our numerical analysis, we first discuss 
the case where /i is purely imaginary, fi = ijl, which is a 
U(l) rotation of the Wilson line: 



V qk (L) 



(30) 



If the overall symmetry were U (3), instead of SU (3), then 
the Vandermonde potential is independent of Jl. For a 
SU(3) theory, however, the Z(3) symmetry only requires 
that the potential is degenerate when Jl — and ±27r/3. 

As Jl represents an ordinary background field, the anti- 
triplet loop is the complex conjugate of the triplet loop. 
Defining 9 as the phase of £, £ = exp(i0)\£\, then £* = 
exp(— i9)\£\, and the Vandermonde potential is a function 
of \£\ and 9. 
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FIG. 2: The expectation values {£} and (I*) as functions of /i 
for h = 0.1 and ~- " 
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When the chemical potential is real, as noted from 1|12|) 
and lfT5)l. {£) and (£*) are unequal but real. Fig. |21 shows 
the expectation values of the loop and of its conjugate 
for a background field ft = 0.1 (again without a gluon 
loop potential). For small ft and fi = 0, an analytical 
discussion of the N = oo potential at the Gross- Witten 
point shows that (£) ~ ft/2, cf. section IIIB in [3^|. From 
Fig. |21 one observes that this remains approximately true 
also for three colors. 

At non-zero \i then, {£) and {£*) split, which is due to 
the imaginary part of the fermion contribution JHJ) to the 
loop action. While (£*) increases monotonically with /i, 
{£) initially decreases from its value at fi = 0, cf. cq. 114fl . 
Finally, both expectation values approach one at large /i, 
in accord with our discussion in Sec. [I] 

In Sec. [H] we saw that in a U(l) model, the in- 
dependence of the expectation values is entirely given 
by a factor exp(±/x), (|19fl . We have checked numeri- 
cally that this is approximately valid for SU(3) when the 
chemical potential is very small. Fig. [21 shows, however, 
that this fails when fi ~ 1. 
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In a matrix model, £3 and % are traces of matrices. 
One could also consider a Polyakov loop model |28j . 
where £3 and % are just scalar fields. To reduce the 
global symmetry from U(l) to Z(3), it is necessary to 
include cubic terms, such as (£ 3 ) 3 + (%) 3 , ©. As for the 
matrix model, one finds that the expectation values of 
£3 and £-g differ when /1 ^ 0. Their exact form depends 
upon the details of the loop potential. 

We now turn to a discussion of the effective potential 
V e ff (£,£*), for real chemical potential. We also include 
the gluon loop potential from eq. (J2J. The solutions of 
the stationarity conditions 



dV eif {£,£*) _ dVeff (£,£*) 



d£ 



d£. 







(31) 



determine the expectation values of the triplet and anti- 
triplet loops, (£3) and (^3}. Due to (|26|l . these equations 
can be rewritten as 



co(C£° 



u*(C£° 



2M 



m z £ 



d£ t 



[V g i (£, 4) +V qk (£,£*)] 



m 2 £° + - e->* 



(32) 
(33) 



These equations have to be solved simultaneously 
with l|24l) . Note that at the stationary point both to and 
lu* are real. Also, these equations, unlike (|T2l and lfT3|) 
above, make it obvious that the expectation values of the 
loops are not determined by cancellations of positive and 
negative contributions: (|32|1 and i|33|l do not involve any 
oscillating functions. 

To show the shape of the effective potential we fix £ — 
£* to its expectation value given in eqs. (|12I13[) or in 
eqs. I|32I33[1 above. We then study V e // as a function of 
the remaining degree of freedom, £ + £„ . 

The behavior of the effective potential with nonzero h 
and /x is customary of a first order transition in a back- 
ground magnetic field. Fig. [3] shows the effective poten- 
tial for h = 1CP 3 and /x = 0, 2, 3, respectively, as a func- 
tion of £ J r£. lf . For each curve, the coupling m 2 — m 2 (/x, h) 
is adjusted to maximize the susceptibility d(£3+£^)/dm 2 . 
For such a weak background field, the first-order phase 
transition persists at fx = 0. As /x increases, the two 
minima of V e // approach each other and the barrier de- 
creases. The first-order phase transition ends in a critical 
point at /x = /xe- The transition is of second order at /xg, 
as the mass of the real part of the triplet loop vanishes. 
From Fig. EI ~ 2.0. As the chemical potential in- 
creases above lie, the mass increases again, and there is 
no phase transition. 

To date, Monte Carlo simulations have been performed 
for large T and small /i 0, 0, |2(J 0] . A matrix model 
predicts that (£) ^ {£*) when fi ^ 0. The first work of 
Allton et al. [n| did not test this directly, but finds that 
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FIG. 3: The effective potential at the peak of the Polyakov 
loop susceptibility for h = 10~ 3 and various fj,. 



{£) changes when the sign of fi is flipped. They did not 
plot (£) or {£*) versus /1, and so did not test the prediction 
that one of these expectation values is not monotonic in 

The analysis of the present paper is most applicable 
for heavy quarks. At fx = 0, the lattice gives us an out- 
line of the phase transition for three degenerate flavors of 
massive quarks. The deconfining transition only persists 



for relatively heavy quarks, m > m 



end 
qk 



Q; 



this 



was estimated to disappear for a pseudo-scalar mass of 
« 1.8 GeV. There is no phase transition for intermediate 
quark masses, with a first order chiral transition appear- 
ing for light quark masses. In all cases at fx = 0, the 
rise in the Polyakov loop appears to coincide with the 
decrease in the chiral order parameter. 

The case of heavy quarks at /i ^ is then similar 
to that of Fig. [3] a first order transition at fi = 0, 
ending in a critical end point at some \xe- (See, also, 
Fig. 1 of 0). For quarks lighter than m^ d , there is 
no deconfining transition, and the correlation length of 
the Polyakov loop decreases monotonically as fj, increases 
from zero. In the plane of fx and T, there may be a crit- 
ical end point at fx ^ 0, where the correlation length for 
the sigma meson diverges ^3 i that for the Polyakov loop 
will remain finite, except from its coupling to the sigma. 

To describe the region of small quark masses, and the 
chiral transition, it is necessary to introduce a chiral or- 
der parameter, and couple that to the Wilson line. A 
mean field approximation can be analyzed in a matrix 
model with two coupled matrices |4l|. Due to the large 
N argument mentioned at the end of Sec. II, it is possible 
that for three colors, the coincidence of the chiral and 
deconfining "transitions", ubiquitous at jx = 0, breaks 
down at some finite value of [i. That is, for intermediate 
quark densities, hadronic matter exists as a Fermi sea of 
"confined", but chirally symmetric, nucleons. 
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IV. LATTICE QCD AND THE FERMION SIGN 
PROBLEM 

We conclude by discussing how the results of the ma- 
trix model may be of use for numerical simulations of 
dense quarks in lattice QCD. 

In Euclidean spacetime, the quark part of the action 

is 

S qk = J d A x i> (p+ ^7° + m) ip . (34) 

We follow the conventions of 39], with p= {da—igAJ)^ 
the covariant derivative for a gluon field A^. In this sec- 
tion, and in contrast to previous notation, here /x is the 
quark chemical potential (not fi/T), and m is the quark 
mass (not m q k). 

We need to use two symmetries. By a combination of 
Hcrmitian conjugation, plus a 75 transformation, 

(det(p+ ^7° + m))* = det{fl- + to) ; (35) 

see, e.g., (13) of . At zero chemical potential, Jp is 
purely anti-Hcrmitian; as Jj) anti-commutes with the ma- 
trix 75, the eigenvalues pair up, and the quark determi- 
nant is real. At nonzero chemical potential, the quark 
operator is a sum of an anti-Hermitian operator, IJ), and 
a Hcrmitian operator, /j-jq. While the eigenvalues form 
pairs with opposite sign, l|33|) shows that the quark de- 
terminant is complex when /j ^ 0. This is the fermion 
sign problem in dense QCD. 

We can perform a charge conjugation transformation 
on the quarks [3^| . This is a change of variables in the 
Grassman integration over the quarks, and so it doesn't 
change the determinant. This gives 

det (lp-\- fxy + m ) = det (fl c - ^7° + m) , (36) 

where lp c = (d^ + igA*^ )7 M is the covariant derivative for 
the charge conjugate gluon field, —A* . By itself, this 
isn't of much help, as we have changed the sign of the 
chemical potential, and turned the gluon field into its 
charge conjugate. In the matrix model, this symmetry is 
manifest in (J5J. 

We now combine these two relations, to obtain: 

det(p c + /l«7 + to) = (dct(^+ M 7 n + ?7i))* . (37) 

This shows that for the same sign of /x, the quark deter- 
minant for charge conjugate gluons is the complex conju- 
gate of the quark determinant in the original gluon field. 
This generalizes what is obvious in the matrix model. 

Thus while the quark determinant in the presence of a 
given gluon field is complex, by adding the contribution 
of the charge conjugate gluons, we immediately obtain a 
partition function whose measure is manifestly real. This 
extends immediately to the lattice. There, gluons live on 
links, with link fields Un = exp(igaA^), where "a" is 
the lattice spacing. A configuration of links is given by 
some set of U^'s; the charge conjugate lattice is simply 



given by replacing each JJ ^ by U*. That one can, in this 
way, obtain a real measure of the functional integral was 
known from the work of the Glasgow group (see, e.g., the 
discussion just before Eq. (8) in the last reference of 0). 

However, all this does is to reduce the problem from 
one of complex phases, to one of real phases. There are 
still configurations in the functional integral with both 
positive and negative weight. This still leaves the prob- 
lem of how to decide whether to sum over configurations 
with both signs. Also, how does one include the effects 
of a Fermi sea of quarks in weighting configurations? 

The matrix model provides clues to both of these ques- 
tions. It is true that configurations of both signs con- 
tribute to the integral of the matrix model. However, at 
zero density, the background Z{N) field which quarks in- 
duce provides an expectation value along a definite direc- 
tion in the complex plane, for real values (this is related 
to the sign of the quark masses). Further, at nonzero 
quark density, the field for the imaginary part of the loop 
has an imaginary coefficient, so that the expectation val- 
ues of both t and I* remain real and positive. We have 
checked that even at nonzero /i, the dominant configura- 
tions of the matrix model are those in which the measure 
is positive. 

It is reasonable to conjecture that this remains true 
with dynamical quarks. This suggests that in Monte 
Carlo simulations, that one accept configurations in 
which the quark determinant is positive, and drop those 
in which it is negative. 

How, then, does one weight by configurations which 
include the effects of a Fermi sea of quarks? This could 
be included by a type of tadpole improvement. Suppose 
that one works out from zero chemical potential, to in- 
creasingly large values. To represent the effects of fi ^ 0, 
one would expand not about the bare link variables, but 
about links equal to the expectation value of the loop. 
For a link going forward, one would use (£); for a link 
going backward, (£*). This will explicitly bias one to 
configurations which include, approximately, the effects 
of the Fermi sea. 

This is supported by numerical simulations of Blum, 
Hetrick, and Toussaint for heavy quarks [lpj]. Using their 
results, de Forcrand and Laliena showed that the 
phase of the (heavy) quark determinant is proportional to 
the phase of the Polyakov loop, times the spatial volume. 

These results illustrate a more general problem. The 
parameters of a matrix model are just numbers. This 
represents, however, a mean field approximation to the 
theory in a spatial volume, V. For example, the back- 
ground field induced by quarks is itself proportional to 
V. For a measure which is always real and positive, this 
is of no concern: even if h ~ V, an error of order one 
is inessential relative to the dominant term, which is 
cxp(— #h) ~ exp(— #V). The integrals which enter at 
nonzero density, however, are those in which the mea- 
sure includes oscillatory terms, as in @, (|12p. and 11311 . 
In this case, it is necessary to determine the phase accu- 
rately not just to ~ V, but to ~ 1! In essence, this is the 
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true fermion sign problem: not that the measure is not 
positive definite, but that one must determine the phase 
of the quark determinant very accurately. We note that 
similar oscillations in the quark determinant have been 
derived, using random matrix theory in the e-regime, by 
Osborn, Splittorff and Verbaarschot |2j| • 

Nevertheless, we suggest that these techniques might 
be of use in numerical simulations of dense QCD on the 
lattice. By their nature, they are most suited for heavy 
quarks, starting from the region of zero density, and 
working out to nonzero density. Even if one accepts con- 
figurations whose overall weight is positive, it is certainly 
necessary to use cluster algorithms to include regions in 
which the phase is negative [Tfil |. 

We have used an effective model which is, implicitly, 
valid only for distances 3> 1/T. When the temperature 
is small, it is also imperative to include fluctuations in 
the expectation values of timelike links, as they wander 



about in (imaginary) time. 

Lastly, these ideas are strongly motivated by the heavy 
quark limit [Tfl Ill1 |. where quarks only propagate up- 
ward in imaginary time. Light quarks also propagate in 
space, so that at nonzero density, one will have to expand 
about modified expectation values for propagation which 
is "forward" or "backward" in proper time. 

There is now a wealth of results available at nonzero 
temperature and small chemical potential [Tsl IT9II20II21I 
E^ . 12^ . This is the first place to test our admittedly 
speculative remarks. 
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